Bypassing the bandwidth theorem with VT symmetry 
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The beat time Tf pt associated with the energy transfer between two coupled oscillators is dictated 
by the bandwidth theorem which sets a lower bound rf pt ~ 1/Suj. We show, both experimentally and 
theoretically, that two coupled active LRC electrical oscillators with parity-time (VT) symmetry, 
bypass the lower bound imposed by the bandwidth theorem, reducing the beat time to zero while 
retaining a real valued spectrum and fixed eigenfrequency difference Slu. Our results foster new 
design strategies which lead to (stable) pseudo-unitary wave evolution, and may allow for ultrafast 
computation, telecommunication, and signal processing. 
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One of the fundamental principles of wave physics is 
the Bandwidth Theorem [1] which in quantum mechanics 
takes the form of the celebrated energy-time Heisenberg 
Uncertainty relation [2]. A direct consequence of this 
principle is the fact that the time for evolution between 
two orthogonal states Tf pt (first passage time) is bounded 
by Tfpt <~ 1/Siu [3, 4]. A basic example where this lower 
bound can be exhibited is the beat time associated with 
the energy transfer between two coupled oscillators. Even 
though the validity of the uncertainty principle is un- 
doubted, it has recently been suggested that possible ex- 
tensions of quantum mechanics invoking non-Hermitian 
VT— symmetric Hamiltonians [5] can generate arbitrar- 
ily fast state evolution referred to as brachistochrone dy- 
namics [6, 7, 9-11]. The main characteristic of this class 
of Hamiltonians % [6, 8-10] is that they commute with 
an anti-linear operator VT, where the time-reversal op- 
erator T is the anti-linear operator of (generalized) com- 
plex conjugation and V is a (generalized) parity operator 
[5] . Examples of such PT-symmetric systems range from 
quantum field theories to solid state physics and classical 
optics [17-34]. Due to the anti-linear nature of the VT 
operator, the cigenstates of % may or may not be eigen- 
states of VT- In the former case, all the eigenvalues of % 
are strictly real and the PT-symmetry is said to be exact. 
Otherwise the symmetry is said to be spontaneously bro- 
ken. In many physical realizations, the transition from 
the exact to the broken VT -symmetric phase is due to 
the presence of various gain/loss mechanisms that are 
controlled by some parameter 7 of % . 

At the same time, the brachistochrone evolution has 
a long standing history and is significant both in theory 
and in application. It was one of the earliest problems 
posed in the calculus of variations, and in the frame- 
work of classical mechanics it dictates "the curve down 
which a particle sliding from rest and accelerated by 
gravity will slip (without friction) from one point to an- 
other in the least time" [35]. The quantum mechani- 
cal brachistochrone problem has recently been revived 
in the emerging fields of quantum computation and sig- 



nal processing where one studies the possibility to use 
dynamical protocols in solving computational problems 
and enhancing signal transport respectively. More specif- 
ically, the quantum brachistochrone problem can be for- 
mulated as follows: Given two orthogonal states and 
I one wants to find the (time-independent) Hamil- 
tonian H (protocol) that performs the transformation 
— > \^ f) = e^ tHt \^i) in the minimal time Tf pt for 
a fixed difference Slj = \(Ef — Ei)/h\ of the eigenval- 
ues Ef,Ei of H [36]. Such a constraint is appropriate 
since a rescaling of the Hamiltonian as H — > AH, with 
A > 1, would make Slu, and hence the transition rates, 
large. This corresponds to the fact that physically only 
a finite amount of resources (e.g. a finite magnetic field, 
bandwidth resolution, etc.) are typically available. The 
quantum brachistochrone was addressed by a number of 
researchers (see for example [3, 4, 36-38]) who show that 
the minimal time to perform the required transformation 
is bounded by the energy-time uncertainty Tf pt <~ j^. 

Here we adopt these ideas into the realm of classi- 
cal wave propagation and engineer a system comprised 
of two active LRC circuits with PT-symmetry having 
an arbitrarily low first passage time Tf pt bypassing the 
lower bound imposed by the Bandwidth theorem. We 
call this limit-breaking wave phenomenon tachistochrone 
passage [39] where the equivalent system of coupled PT- 
symmetric oscillators with the same bandwidth Slu can 
be obtained if one oscillator experiences attenuation with 
rate 7 while its partner experiences an equivalent amplifi- 
cation rate 7. Depending on the application point of the 
incident excitation, we observe unidirectional accelerated 
signal/energy transport, where the time Tf pt <~ I/7 can, 
in principle, become arbitrarily short. Although these 
results rely on the fact that the generated dynamics is 
non-unitary, we mark that due to the PT-symmetry, the 
corresponding eigenfrcquencies to can be real, thus guar- 
anteeing the stability of our system [40]. Our results 
foster new design strategies based on active elements [41] 
with VT symmetric arrangements that may allow for ul- 
trafast computation, telecommunication, and signal pro- 



FIG. 1: Electronic implementation of a 'PT-symmetric dimer. 
The coils are inductively coupled with Vi and V2 providing ac- 
cess to the system variables. The switch S asserts the initial 
condition. A resistor provides the damping (green side) while 
a negative resistance gain element (red side) is implemented 
by feedback from an LM356 voltage doubling amplifier. Cur- 
rent flows in the direction of the arrows proportional to the 
respective voltages Vi and V%. 



cessing. 

Our system, shown in Fig. 1, consists of a pair of cou- 
pled LC circuits, one with amplification and the other 
with equivalent attenuation. The circuit was shown in 
Ref. [31] to be a simple realization of a PT-symmetric 
dimer. Each inductor is wound with 75 turns of #28 cop- 
per wire on 15cm diameter PVC forms in a 6 x 6mm loose 
bundle for an inductance of L = 2.32 mH. The coils are 
mounted coaxially with a bundle separation adjusted for 
the desired mutual inductance M. The isolated natural 
frequency of each coil is ljq = 1/y/LC = 2 x 10 5 s _1 . 

The actual experimental circuit deviates from Fig. 1 in 
the following ways: (1) A resistive component associated 
with coil wire dissipation is compensated by an equiv- 
alent gain component applied to each coil; (2) A small 
capacitance trim is included to aid in circuit balancing; 
and (3) Additional LM356 voltage followers are used to 
buffer the voltages V\ and V 2 , captured with a Tektronix 
DPO2014 oscilloscope. 

The linear nature of the system requires a balance of 
VT symmetry only to the extent that component stabil- 
ity over time allows for a measurement. All circuit modes 
either exponentially grow to the nonlincarity limit of the 
buffers, or shrink to zero. Transient data is obtained 
respecting these time scales. 

Kirchhoff's laws lead to the following set of equations 
for the charge Q\ (Q 2 ) on the capacitor corresponding to 
the amplified (lossy) side: 



d 2 Q 1 
dr 2 

d 2 Q 2 
dr 2 



—aQi + fiaQ 2 + 7 
HaQi - aQ 2 - 7 



dQi 
dr 
dQ 2 
dr 



(1) 



where r = to t, a = 1/(1 — ^i 2 ) > 1, 7 = Rr X yjL/C is the 
gain/loss parameter, and fi = M/L is the rescaled mutual 
inductance. Inspection of Eqs. (1) reveals that they are 
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FIG. 2: (Color online) Parametric evolution of the "PT-dimer 
eigenfrequencies vs. the gain/loss parameter 7, for different 
mutual inductances [i (indicated in the figure) . The solid cir- 
cles are experimental data. The borders of the highlighted 
area (marked by a black dashed line) indicate a constant 
bandwidth path with fixed frequency difference Suj/uiq = 0.36. 



invariant under a combined parity (i.e. Qi Q2) and 
time-reversal (i.e. t — > —t) transformation. 

The theoretical analysis (see reference [42]) relies on a 
Liouvillian formulation of Eqs. (1) which take the form 
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where * = (Qi,Q 2 , Qi,Q 2 ) T - Eq. (2) can be interpreted 
as a Schrodinger equation with non-Hermitian effective 
Hamiltonian H c ff = iC This Hamiltonian is symmet- 
ric with respect to generalized VqTq transformations, i.e. 
[PqTo, H eS ] = 0, where 
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and (j x is the Pauli matrix, 1 is the 2x2 identity matrix, 
and K, denotes the operation of complex conjugation. By 
a similarity transformation 1Z, 



K = 



j b + c b + c i — 
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H c g can be related to a transposition symmetric, 
VT— symmetric Hamiltonian H. Specifically, 



H = H T = KH eff K-\ 
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FIG. 3: (Color online) A representative tachistochrone passage wave propagation for a circuit with a bandwidth constraint 
Soj/loo = 0.36. (a) A typical temporal dynamics of the displacement current -Zi(t) when the initial condition corresponds to an 
excitation of the circuit n = 1 (gain side); (b) the same for current 72 (t) when the initial condition is an excitation at circuit 
n — 2 (loss side). The black lines are numerical simulations for the passive 7 = circuit, while the red and green lines are 
numerical simulations for the active circuit with 7 = 0.24 and /1 = 0.39. The measurements are indicated with circles. The 
first passage time rf pt is indicated with a cyan (black) arrow for the active (passive) dimer. (c) The extracted Tfyt versus the 
gain/loss parameter 7. The green circles correspond to an initial condition 72(0) = 1 (loss side) while the red circles to an 
initial condition 7i(0) = 1 (gain side). The solid lines indicate the theoretical result Eq. (6). The black dashed line indicates 
the beat time for a passive (7 = 0) system, (d) The first passage time Tfp t versus the gain/loss parameter 7 as it is measured 
from the energy exchange condition (see text). Open circles are the experimental data while lines (of similar color) correspond 
to the theoretical prediction Tfpt = Tf p t/2. The black dashed line denotes r f p t for 7 = 0. 



where 
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and b = v / (a + a 1 / 2 )/2, c = -^{a - a 1 / 2 )/2. This al- 
lows us to make contact with the brachistochrone studies 
of Refs. [6, 8, 9, 15]. 

The eigenfrequencies uj\^ of system (1) (or equiva- 
lently Eq. (2)) are shown as functions of the gain/loss 
parameter 7, in Fig. 2 (solid lines). For 7 < 7737- = 
l/\/l — fJ- — l/VT+/i the system is in the exact phase 
and thus the eigenfrequencies are real [31]. We investi- 
gate the signal/energy tachistochrone passage under the 
constraint of fixed bandwidth Soj — oj\ — 0J2 ■ Experimen- 
tally, the Soj constraint is implemented in the LRC dimer 
through adjustment of the mutual inductance. 

Figure 2 shows both theoretical and experimental re- 
sults for the parametric evolution of the eigenfrequencies 
in the exact phase, for various [i values. The black dashed 
lines in Fig. 2 illustrate a path for fixed SojjojQ = 0.36 
through the family of eigenfrequencies associated with 
different mutual inductances \x. 

Eqn. (2) can be solved either analytically or via di- 
rect numerical integration in order to obtain the tem- 
poral behavior of the capacitor charge Q n {j) an d the 
displacement current 7„ (r) in each of the two circuits of 
the 'PT-symmetric dimer. For the investigation of the 
tachistochrone wave evolution, we consider an initial dis- 
placement current in one of the circuits with all other 
dynamical variables zero. The first passage time Tf pt is 
then defined as the time interval needed to reach an or- 
thogonal state. In our experiments this corresponds to 
the condition that the envelope function of the current at 
the initially excited circuit is zero. We find that the first 
passage time is asymmetric with respect to the initially 



excited circuit. Specifically we have that 
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(see reference [42]) where the + sign corresponds to an 
initial condition starting from the gain side while the — 
sign corresponds to an initial condition starting from the 
lossy side. For 7 ^> Soj, Eq. (6) takes the limiting values 
Tfpt ~ 2-k/Soj and Tf pt f=s 2/7 respectively. The latter case 
indicates the possibility of transforming an initial state to 
an orthogonal final one, or in more practical terms, trans- 
ferring energy from one side to the other, in an arbitrarily 
short time interval. In the opposite limit of 7 = 0, we re- 
cover for both initial conditions the Anandan-Aharonov 
lower bound for the first passage time Tf pt = it/Soj [3]. 
This is the time for which energy is transferred from the 
initial circuit to its partner according to the constraint 
of the Bandwidth theorem. 

Geometrically, one can understand the relation (6) in 
the following way: the time required for the evolution 
between two states induced by a Hermitian Hamiltonian 
is proportional to the length of the shortest geodesic 
connecting the two states in projective Hilbert space 
[3]. Non-Hermitian VT- symmetric Hamiltonians in the 
exact PT-symmetric domain can be similarity mapped 
to equivalent Hermitian Hamiltonians. Under such a 
similarity mapping the corresponding projective Hilbert 
space undergoes a deformation obtaining a nontrivial 
metric. This results in an effective contraction or di- 
lation of the corresponding geodesic and with it of the 
corresponding evolution time [8, 9, 16]. 

In Fig. 3 we present some typical measurements for 
the temporal behavior of displacement currents. In sub- 
figure 3a we show |7i(t)| for an initial condition corre- 
sponding to the case 7"i(0) = 1 with all other dynamical 
variables zero. The case where the initial current exci- 
tation is at the lossy side i.e. 72(0) = 1, is shown for 
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contrast in Fig. 3b. In both cases, agreement between 
the experiment (circles) and the simulations (lines) is ob- 
served. For comparison, we also report with black line 
the temporal behavior of the displacement current for the 
case of a passive circuit (i.e. 7 = 0) with the same 5u>- 
constraint. We observe that the orthogonal target state 
is reached faster (or slower) depending on whether the 
initial excitation is applied to the lossy (or gain) side. 

The above results can be verified in more cases by 
changing the inductive coupling /x and gain/loss param- 
eter 7, while keeping constant the frequency difference 
Slo = lu 2 — u>\ . A summary of our measured Tf p t versus 
7 is presented in Fig. 3c. The experimental data show 
agreement with the theoretical prediction Eq. (6). 

A parallel analysis pertains directly to the study of 
the energy transport from one side to another. Using 
the same initial conditions as above we investigate the 
temporal behavior of the energies 



E n (r) 



lQl + 1 LI* 



(7) 



of each n = 1,2 circuit. The first passage time can be 
defined as the time for which the two energies become 
equal for the first time i.e. Ei{t^ x ) — E 2 (r^ t ). For pas- 
sive (i.e. 7 = 0) coupled circuitry, this time is half of 
the beating time r f ^ t = Tf pt (7 = 0)/2 and it is insensi- 
tive to the initial preparation. In contrast, for the active 
PT-symmetric dimer of Fig. 1, we find that the energy 
transfer from the lossy (gain) side to the gain (lossy) one, 
is faster (slower) than the corresponding passive system 
with the same 6u. In Fig. 3d, we summarize our mea- 
surements for the r f p t versus 7 under the constraint of 
fixed frequency bandwidth 5u. A similar behavior as the 
one found for the displacement current is evident. 

Our results open a new direction towards investigating 
novel phenomena and functionalities of PT-symmetric 
arrangements with pseudo-unitary spatio-temporal evo- 
lution. Along these lines, we envision PT-symmetric 
(nano)-antenna configurations and metamaterial or op- 
tical microresonator arrays with unidirectional ultra-fast 
communication capabilities. These structures also have 
potential applications as delay lines, buffers and switches. 
Questions like the effects of non-linearity or the topolog- 
ical complexity of the PT-symmetric structures in the 
tachistochronc dynamics are open and offer new exciting 
opportunities, yet to be discovered. 
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SUPPLEMENTAL MATERIAL 

From the eigenvalue equation {H e ff — Wfc)Sfc = with 
k = 1, 2, 3, 4 and the corresponding characteristic polyno- 
mial io\ — (2a — j 2 )ujl + a = one obtains 07^ = ±\/Tt. 
w 2 ,3 = where 
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a — UJi 
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ajj, 



iLo^e 



a k e 
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For gain/loss parameters 7 e [0, j-pr] the eigenvalues are 
purely real, u>/. € K; it holds 4>k £ K so that PoT^fe = 
and the PoT — symmetry is exact. For 7 > "f-pT the 
eigenvalues u>k are not real, but pairwise complex con- 
jugate; one finds (f>k ^ R so that PoTSfe 9^ Sfc and 
the PoT - symmetry is spontaneously broken with PoT 
phase transition at 7 = 7737-. The Hamiltonian H c ^ can 
be brought into the more symmetric form Eq. (5) via 
a similarity transformation 1Z given by equation (4). To 
describe the dynamics in the sector of exact PoT symme- 
try, 7 € [0, 7pt]i where u>i = —lo^ 1 lj 2 = —1^3, W/t € K we 
start from an ansatz ^(t) = X)fe=i e ~ %UkT -^-k^-ki Ak E C 
and impose the experimentally required reality constraint 
*(t) e K 4 as AiSx = A4S4, A 2 E 2 = A 3 E 3 . The solu- 
tions of the evolution equation (1), (2) with initial con- 
dition on the gain side ty(r = 0) = (0, 0, 1, 0) T are 
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The first passage time Tf pt required for the evolution 
from the initial state *(r = 0) = (0,0,1,0) T to an 
orthogonal final state ^(r = Tf pt ) = (ipi, ip 2 , 0, tp4) T , 
(5'(0)| v I , (Tf pt )) = with regard to the slowly evolv- 
ing enveloping amplitude follows from tp3(r) = 

V /l + 3&\A+5 Sin [ SuJT+ t- 52 ] sm [^±^] as 



Tf p t = (5 2 - 5\)/5ui. 



(A-3) 



Due to the invariance of equations (1) under simulta- 
neous action of Qi <^ Q 2 and 7 — > —7 the evolu- 
tion from the lossy side *(t = 0) = (0,0,0, 1) T to 
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^(Tfpt) = (^l? ^2; "03i 0) T requires a Tf p t obtainable from 
(A-2) by sign change 7 — > —7. Explicitly this yields Eqs. 
(6). 
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